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According to the Banks-Casher formula the chiral order parameter is directly 
related to the spectrum of the Dirac operator. In this lecture, we will argue 
that some properties of the Dirac spectrum are universal and can be obtained 
from a random matrix theory with the global symmetries of the QCD partition 
function. In particular, this is true for the spectrum near zero on the scale 
of a typical level spacing. Alternatively, the chiral order parameter can be 
characterized by the zeros of the partition function. We will analyze such 
zeros for a random matrix model at nonzero chemical potential. 

1. Introduction 

Many phenomena in nuclear physics, as for example the lightness of the 
pion mass and the absence of parity doublets, can be explained by the as- 
sumption that chiral symmetry is spontaneously broken. This assumption has 
been confirmed by numerous lattice QCD simulations (for a review see 0,0]). 
However, these studies also show that chiral symmetry || is restored at a 
critical temperature of T c « 140 MeV. In spite of steady progress Q, the 
situation at nonzero baryon number density is much less clear ||. It seems 
that the quenched approximation does not work |5|J|, and the phase of the 
fermion determinant makes unquenched simulations virtually impossible. 

The order parameter of the chiral phase transition is the chiral condensate. 
It is directly related to the spectral density of the Euclidean Dirac operator [[7] . 
One of the questions we wish to address is to what extent the Dirac spectrum 
shows universal features which can be obtained from a Random Matrix Theory 
(RMT) with the global symmetries of the QCD partition function (chiral Ran- 
dom Matrix Theory (chRMT)). As is well-known from the study of complex 
systems || , only correlations on a scale set by the eigenvalue spacing are given 
by RMT. Such correlations may be important in mesoscopic systems. Typical 
examples are a finite nucleus ||, small metallic particles |TIJ, quantum dots 



and disordered wires (see 0] f° r a review). In particular, universal conduc- 



tance fluctuations have been understood in the framework of RMT 111211 . In 
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lattice QCD simulations, with a mesoscopic number of degrees of freedom, we 
expect to observe similar phenomena. In particular, in the mesoscopic range 
of QCD, for box size L given by JTS] 1/A <C L <C we expect to obtain 

exact results from RMT (A is a typical hadronic scale and m n is the pion 
mass) . 

In spite of its success in explaining mesoscopic phenomena, I wish to em- 
phasize that there is no universality on a macroscopic scale and that RMT 
cannot be used to obtain quantitative predictions at this scale. (Note, how- 
ever, [p^ -[l7|.). This does not imply that RMT cannot be useful for the study 



of macroscopic phenomena. They have widely been used as schematic models 
for disorder, e.g., Anderson localization and the Gross- Witten model ||19| . 



Because of problems in simulating QCD at a nonzero baryon number density, 
chRMT is an ideal laboratory to address this problem. One important success 
is the understanding of the nature of the quenched approximation ||. Below, 
we will discuss the phase structure of chRMT by means of the distribution of 
Yang- Lee zeros |20| . 



2. The Chiral order parameter 

The order parameter of the chiral phase transition, (ipip)., is nonzero only 
below the critical temperature. As was shown in |7j (ipip) is directly related to 
the eigenvalue density of the QCD Dirac operator per unit four-volume 

E.|(«| = ^f» (1) 

It is elementary to derive this relation. The Euclidean Dirac operator for gauge 
field configuration is given by D = 7^(9^ + iA^). For Hermitean gamma 
matrices D is anti-hermitean with purely imaginary eigenvalues, D<pk = i\k<pk, 
and spectral density given by p(A) = X^<KA ~ Because {75,/}} = 0, 
nonzero eigenvalues occur in pairs ±Afc. In terms of the eigenvalues of D the 
QCD partition function for Nf flavors of mass m can then be written as 



([[(% + ™ 2 ) Nf ), (2) 



where the average (•) is over all gauge field configurations weighted according 
to the Euclidean action. 

The chiral condensate follows immediately from the partition function (|2|), 

m = J_ Smlog z (ro) = ^ ? _J^_). (3 ) 
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If we express the sum as an integral over the spectral density, and take the 
thermodynamic limit before the chiral limit so that we have many eigenvalues 
less than m we recover (|I|) (Notice the order of the limits.). 

Another way to characterize the chiral condensate is via the zeros of the 



partition function |2l|]22|| . For a finite number of degrees of freedom Z(m) can 



be factorized as Z{m) ~ Y\k( m ~ m k)- The chiral condensate is then given by 

(#) = T^rd m log Z(m) = i £ (4) 
VJSf V ~r m — rrik 

In the chirally broken phase zeros are located on a segment of the imaginary 
axis that includes m = 0. In the thermodynamic limit they coalesce into a 
cut and the chiral condensate shows a discontinuity each time m crosses this 
cut. In the chirally symmetric phase, we expect to find a cut away from the 
imaginary axis. 

The eigenvalues of the Dirac operator D = 7^(9^ + iA^) + 7oyU at nonzero 
chemical potential are scattered in the complex plane (see for example ||). 
For a finite number of degrees of freedom the partition function, Z(m,fi), is 
a polynomial in m and fi. The condensate can be related to either to the 
spectral density or to the zeros of Z(m,/j,). An alternative order parameter 
is the baryon density n B = log Z(m, n)/NfV. We can differentiate with 
respect to fi before or after averaging over the gauge fields. In the first case 
the baryon number density follows from the spectral density of 70 (-D + m) with 
eigenvalues scattered in the complex plane. In the second case the baryon 
number density follows from the zeros of Z(m,fi) in the complex fi plane. 

3. The Dirac Spectrum 

An important consequence of the Bank-Casher formula ([I]) is that the eigen- 
values near zero virtuality are spaced as AA = l/p(0) = n/J^V. This should be 
contrasted with the eigenvalue spectrum of the non-interacting Dirac operator. 
Then p free (A) ~ VX 3 which leads to an eigenvalue spacing of AA ~ 1/V 1 ' 4 . 
Clearly, the presence of gauge fields lead to a strong modification of the spec- 
trum near zero virtuality. Strong interactions result in the coupling of many 
degrees of freedom leading to extended states and correlated eigenvalues. On 
the other hand, for uncorrelated eigenvalues, the eigenvalue distribution fac- 
torizes and we have p(A) ~ \ 2N f +1 , i.e. no breaking of chiral symmetry. 

Numerous studies have shown that spectral correlations of complex systems 
on a scale set by the level spacing are universal, i.e. they do not depend on the 
dynamics of the system and are completely determined by symmetries. Be- 
cause the QCD Dirac spectrum is symmetric about zero, we have two different 
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types of eigenvalue correlations: correlations in the bulk of the spectrum and 
spectral correlations near zero virtuality. In the context of chiral symmetry we 
wish to study the spectral density near zero virtuality. Because the eigenvalues 
are spaced as 1/HV it is natural to introduce the microscopic spectral density 

ps ^ = ^^k p (vf)- (5) 

The dependence on the macroscopic variable S has been eliminated and there- 
fore ps(u) is a perfect candidate for a universal function. 

4. Spectral universality 

Spectra for a wide range of complex quantum systems have been studied 
both experimentally p3|-|25H and numerically P,p6,p7||. One basic observation 



has been that the scale of variations of the average spectral density and the 
scale of the spectral fluctuations separate. This allows us to unfold the spec- 
trum, i.e. we rescale the spectrum in units of the local average level spacing. 
The fluctuations of the unfolded spectrum can be measured by suitable statis- 
tics. We will consider the nearest neighbor spacing distribution, P(S), the 
number variance, S 2 (n), and the A 3 (n) statistic. The number variance is de- 
fined as the variance of the number of levels in a stretch of the spectrum that 
contains n levels on average, and A 3 (n) is obtained from £ 2 (n) by averaging 
over a smoothening kernel. 

These statistics can be obtained analytically for the invariant random ma- 
trix ensembles (see |8|j28[|) defined as ensembles of Hermitean matrices with 
independently distributed Gaussian matrix elements. Depending on the anti- 
unitary symmetry, the matrix elements are real, complex or quaternion real. 
The corresponding Dyson index is given by f3 — 1, 2, and 4, respectively. The 
nearest neighbor spacing distribution is given by P(S) ~ S 13 exp(— apS 2 ). The 
asymptotic behavior of X 2 (n) and A 3 (n) is given by £ 2 (n) ~ (2/ir 2 /3) log(n) 
and A 3 (n) ~ /3£ 2 (n)/2. For uncorrelated eigenvalues one finds that P(S) = 
exp(— 5*), £ 2 (n) = n an d A 3 (n) = ra/15. Characteristic features of random ma- 
trix correlations are level repulsion at short distances and a strong suppression 
of fluctuations at large distances. 

Numerous studies have shown that the spectral correlations of a classically 
chaotic systems are given by RMT. This conjecture has been strengthened by 
recent analytical arguments p9[- j31"|1 and universality arguments [|32" |. 



5. Chiral random matrix theory 

In this section we will introduce an instanton liquid inspired |33j RMT for the 



QCD partition function. In the spirit of the invariant random matrix ensembles 
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we construct a model for the Dirac operator with the global symmetries of 
the QCD partition function as input, but otherwise Gaussian random matrix 
elements. The chRMT that obeys these conditions is defined by P^-^B 

Z; = I DW I] det(X> + m/Je-^^, with V = 
f=i 

and W is a n x m matrix with v — \n — m\ and N = n + m. The matrix 
elements of W are either real ((3 = 1, chiral Gaussian Orthogonal Ensemble 
(chGOE)), complex (f3 = 2, chiral Gaussian Unitary Ensemble (chGUE)), or 
quaternion real ((3 = 4, chiral Gaussian Symplectic Ensemble (chGSE)). 

This model reproduces the following symmetries of the QCD partition func- 
tion: i) The Ua,(1) symmetry. All nonzero eigenvalues of the random matrix 
Dirac operator occur in pairs ±A. ii) The topological structure of the QCD 
partition function. The Dirac matrix has exactly \u\ = |n — m\ zero eigenval- 
ues. This identifies v as the topological sector of the model. Hi) The flavor 
symmetry is the same as in QCD ||37|| . iv) The chiral symmetry is broken 
spontaneously with chiral condensate given by £ = limA^oo np(0)/N. (N is 
interpreted as the (dimensionless) volume of space time.) v) The anti- unitary 
symmetries. For fundamental fermions the matrix elements of the Dirac oper- 
ator are complex for N c > 3 (/3 = 2) but can be chosen real for N c = 2 (j3 = 1) . 
For adjoint fermions they can be arranged into real quaternions (/3 = 4). 

The ensemble of matrices in (|6]) is also known as the Laguerre ensemble. 
Note that its spectral correlations in the bulk of the spectrum are given by the 
invariant random matrix ensemble with the same value of (3 ||38|| . Both types 



of microscopic correlations are stable against deformations of the ensemble. 
This has been shown by a variety of different arguments [|39|-|42" . 



Below we will discuss the microscopic spectral density. For N c = 3, Nf 
flavors and topological charge v it is given by 



11 / \ 

Ps{u) = 2 [ J a( u ) ~ Ja+l{u)Ja-l{u)) , (7) 

where a = Nf + v. The result for N c = 2, which is more complicated, is given 
43| , and the result for the symplectic ensemble is derived in [H3 . 



m 



Together with the invariant random matrix ensembles, the chiral ensembles 
are part of a larger classification scheme. As pointed out in [IJJ, there is a one 
to one correspondence between random matrix theories and symmetric spaces. 

6. Lattice QCD results 

Recently, Kalkreuter [46| calculated all eigenvalues of the lattice Dirac op- 



erator both for Kogut-Susskind (KS) fermions and Wilson fermions for lattices 
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as large as 12 4 . In the the case of 577(2) the anti-unitary symmetry of the KS 
and the Wilson Dirac operator is different [37,I8|. For KS fermions the Dirac 
matrix can be arranged into real quaternions, whereas the Hermitean Wilson 
Dirac matrix <-y 5 _D Wllson can be chosen real. Therefore, we expect that the 
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eigenvalue correlations are described by the GSE and the GOE, respectively 
48(1 - In Fig. 1 we show results for 0,2(11), A^(n) and P(S). The results for 
KS fermions are for 4 dynamical flavors with ma = 0.05 on a 12 4 lattice. The 
results for Wilson fermion were obtained for two dynamical flavors on a 8 3 x 12 
lattice. Other statistics are discussed in f[49 




Lattice studies of the microscopic spectral density are in progress and pre- 
liminary results are promising IfjOfl . However, an alternative way to probe 
the Dirac spectrum was introduced by the Columbia group ||51|| . They stud- 



ied the valence quark mass dependence of the Dirac operator, i.e. S(m) = 
jj J d\p(X)2m/ (A 2 + m 2 ), for a fixed sea quark mass. In the mesoscopic range, 
the valence quark mass dependence can be obtained analytically from the mi- 
croscopic spectral density (0) J52| , 

3^ = x(I a (x)K a (x) + I a+1 (x)K a ^(x)), (8) 



where x = mVYi is the rescaled mass and a = Nf + v. In Fig. 2 we plot this 
ratio as a function of x for lattice data of two dynamical flavors with mass 
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ma = 0.01 and N c = 3 on a 16 3 x 4 lattice. We observe that the lattice data 
for different values of (3 fall on a single curve. Moreover, in the mesoscopic 
range this curve coincides with the random matrix prediction for Nf = v = 0. 
Apparently, the zero modes are completely mixed with the much larger number 
of nonzero modes. For eigenvalues much smaller than the sea quark mass, we 
expect to see the Nf = eigenvalue correlations. 

7. Chiral random matrix model at nonzero chemical potential 

At nonzero temperature and chemical potential the random matrix Dirac 
operator in |G] is given by |)3"| -|5I)| 

r>-( tW + in T + fi\ 

[ iW* + iQ T + H ) ' K > 

where — T ® n (2n + 

Inspired by p6[, the simplest model is obtained by keeping only the lowest 



Matsubara frequency |[53|j54]| . We wish to stress that this model is a schematic 
model of the QCD partition function. Below, we will discuss a model with 
Qt absorbed by the random matrix and /i / 0. Then the eigenvalues of 
T> are scattered in the complex plane. In the quenched approximation its 
distribution was obtained analytically from the Nf limit of a partition 
function with the determinant replaced by its absolute value. To this end the 
RMT partition function was rewritten in terms a cx-model amenable to a saddle 
point approximation. The cx-model shows a second order phase transition at 
the boundary of the spectrum leading to a vanishing curvature and a diverging 
two point function. This was confirmed by an explicit calculation of this two- 
point function in [fjTj], 

In the remainder of this section we consider the unquenched partition func- 
tion for one flavor. Using a multi-precision package [|58[| , we have calculated 
the (Yang-Lee) zeros of the partition function in the complex fi and m plane 
for values of n as large as 192. Results for n = 192 are shown in Fig. 3. Notice 
that the zeros fall on a curve |5!| . 

From a saddle-point analysis it can be shown that for zero mass the model 
shows a first order phase transition along the curve Re[yU 2 +log(yU 2 )] = — 1 P,|5Ti|. 
The discontinuity of ub requires that the zeros of the partition function fall 
along this curve (see left upper figure). At the endpoints (stars) two different 
solutions of the saddle point equation (a cubic equation [Q) coalesce. All 
other curves can be obtained from a saddle-point analysis as well. A schematic 



picture of the phase structure in the complex m plane is shown in |61 . 
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Fig. 3. Zeros of the partition function in the complex m and ji plane. 



8. Conclusions 

We have shown that the microscopic correlations of the QCD Dirac spectrum 
can be explained by RMT and have obtained an analytical understanding of the 
distribution of the eigenvalues closest to zero. We have given an extension of 
this model to nonzero temperature and chemical potential. Its phase structure 
has been mapped out unambiguously by means of Yang-Lee zeros. 
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